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ABSTRACT 

We have used time-resolved spectroscopy to measure the colour dependence of pul- 
sation amplitudes in the DAV white dwarf G 29-38. Model atmospheres predict that 
mode amplitudes should change with wavelength in a manner that depends on the 
spherical harmonic degree £ of the mode. This dependence arises from the convolution 
of mode geometry with wavelength-dependent limb darkening. Our analysis of the six 
largest normal modes detected in Keck observations of G 29-38 reveals one mode with 
a colour dependence different from the other five, permitting us to identify the £ value 
of all six modes and to test the model predictions. The Keck observations also show 
pulsation amplitudes that are unexpectedly asymmetric within absorption lines. We 
show that these asymmetries arise from surface motions associated with the non-radial 
pulsations (which are discussed in detail in a companion paper). By incorporating sur- 
face velocity fields into line profile calculations, we are able to produce models that 
more closely resemble the observations. 
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Despite observations over a broad range of wavelengths 
and via numerous techniques, the star G 29-38 remains an 
enigma. It is the third brightest ZZ Ceti known (V = 13.05), 
and has pulsation amplitudes among the largest measured 
for these variables (up to 6% modulations at optical wave- 
lengths) . It is the most extensively observed large amplitude 
ZZ Ceti star, having been t he su bject of two global obse rving 
campaigns (Winget et al. 1990; Kleinman et al. 1994) and 



G 29-38, and perhaps for all the large amplitude ZZ Ceti 
stars, because measuring mode periods is a prerequisite for 
measuring mass and internal structure using stellar seismol- 
ogy. The only remaining obstacle is to identify the spherical 
harmonic degree £ and radial order n of the modes detected, 
so that their periods can be compared to those of like eigen- 
modes in structural models of white dwarf stars. 



of cou ntless single-site time series measurements (Kleinman 
et al. 1998 ). It has also been the target of several infrared 
and radial velocity studies following th e det ection of an in- 
frared excess by Zuckerman & Becklin ( 1987 ). Nevertheless, 
we still have neither an unambiguous asteroseismological so- 
lution for the star, nor an explanation for the source of the 
mysterious dust apparently responsible for its excess emis- 
sion in the infrared. 



The pattern identified by Kleinman (1995) is sufficiently 
rich that he was able to attempt mode identification using 
the sa me t echniques applied successfully to DO V (W inget 



Recently, Kleinman ( 1995| , also Kleinman et al. 



1995 



et al. p.991] ) and DBV pulsators (Winget et al. |1994J). He 
searched for the (roughly) equal period spacings and fre- 
quency splittings that signify rotationally split non-radial g- 
modes. His attempt was a measured success; he found that 
the pattern of modes was sensible if interpreted as a sequence 
of mostly £ = 1 modes. Unfortunately, his analysis could not 
assure that any individual mode was £ = 1, nor did it allow 
an unambig uous c omparison to structural models (Bradley 
& Kleinman Il996h. 



analyzed all of the optical time-series photometry of this 
object, and showed that in spite of the changing character 
of the pulsation spectrum each season, there is a stable set 
of recurring modes. This is an important breakthrough for 



With this in mind, and having available to us the bright 
portion of a night allocated to very faint sources at the Keck 
II telescope, we decided to test a method for identifying the 
degree of pulsation modes using G 29-38 as our subject. The 
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method was inspired by the work of Robinson et al. ( [1995 ), 
who used Hubble Space Telescope high speed photometry 
in the ultraviolet to measure £ for modes in the star G 117- 
B15A. Their technique exploited the sensitivity of ZZ Ceti 
mode amplitudes to the wavelength of the observations. At 
all wavelengths, the observed amplitudes are diminished by 
cancellation between surface regions with opposite pulsation 
phase. In the ultraviolet, stronger limb darkening changes 
the character of this cancellation and observed mode am- 
plitudes differ from their optical values in a manner that 
depen ds on i. The models calculated by Robinson et al. 
(1995) to explore differences between ultraviolet and opti- 
cal also show amplitude changes within optical absorption 
lines. The character of these changes likewise depends upon 
£, providing the potential to determine £ from optical spec- 
troscopy. 

In an attempt to measure the line profile variations of 
G 29-38 and use them for £ identification, we acquired over 
four hours of time-resolved spectroscopy using the Keck I I 
Low Resolution Imaging Spectrometer (Oke et al. 1995). 
We have described these obser vations in a companion pa- 
per (Van Kerkwijk et al. 1999[ Paper I), and presented an 
analysis of the periodicities present in the total flux and 
line-of-sight velocity curves. Prior to our Keck observations, 
the velocity variations associated with ZZ Ceti pulsations 
had never been detected. Their presence significantly com- 
plicates the models required to fully understand our data, 
but also increases the amount of valuable information we 
can hope to extract. 

In this paper we present our analysis of the line pro- 
file variations of G 29-38. We begin in with an analysis 
of the average spectrum, which has a signal to noise ratio 
higher than usual for white dwarf spectra. In §^| we present 
the amplitudes and phases as a function of wavelength for 
the largest modes and compare them to each other and to 
t heore tical models like those calculated by Robinson et al. 
( 1995 ). On this basis alone it will be clear that we can un- 
ambiguously identify I for these modes. It will also be clear 
that good quantitative fits to the data will require improve- 
ments to model atmospheres. Nonetheless, we will continue 
in §^ by incorporating the velocity field associated with the 
pulsations into the models. The improved models help to 
constrain other pulsation properties, such as the velocity 
amplitude of motions at the stellar surface. We present our 
conclusions in ffcl 



2 THE MEAN SPECTRUM 

By averaging together all of our time-series spectra, we 
have constructed a mean spectrum with extremely high sig- 
nal to noise ratio. This spectrum shows Balmer lines (H/3 
through Hi), the Ca II A3933 resonance line, and a hint of 
Mg II A4481 (Paper I). Metal lines were fir st dis covered in 
G 29-38 by Koester, Provencal, & Shipman (1997), who also 



found iron lines in the ultraviolet spectrum. The presence 
of me tals in a DA spectrum is unusual, and Koester et al. 
( 1997 ) attribute them to the accretion of dust, supporting 
the notion that the infrared excess in the spectrum of G 29- 
38 is caused by reprocessing of light by circumstellar dust 
grains. 

We have fitted the Balmer lines in our mean spectrum 
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Figure 1. The best overall fit to the line profiles H/3 through Hll 
in the average spectrum of G 29-38. The model (dotted line) has 
log g = 8.05 and T cff = 11, 850 K. The dip on the blue side of He 
is the Ca II A3933 line. 



using a grid of model spectra kindly provided by D. Koeste r 
(for a recent description, see Finley, Koester, & Basri 1997). 
The models consist of tabulated intensities, 7^, at 9 limb 
angles, /i = cos(6), for a grid of atmospheres with effective 



temperatures spanning the ZZ Ceti instability strip and with 
gravities from log g — 7.50 to 8.75. The atmospheres were all 
calculated using the ML2, a = 0.6 prescription for convec- 
tion, which yields consistent fits for ZZ Ceti star s ove r the 
broadest range of wavelengths (Bergeron et al. 1995). To 



compare the models to our spectrum we integrated 7a over 
the visible hemisphere and compared the resulting Balmer 
lines to those we observed by normalising both model and 
data to fixed continuum points. This m ethod is similar to 
that used by Bergeron, Saffer & Liebert (1992), but le ss so- 
phisticated than the procedure Bergeron et al. ( |1995| ) used 
for their analysis of ZZ Ceti spectra. 

Figure [l] shows the observed Balmer lines along with the 
best fitting model, which has T c g = 11,850 K and loggr = 
8.05. These valu es are close to those publi shed f or G 29-38 by 
Bergeron et al. (1995) and Koester et al. (1997). The former 
found T aS = 11,820 K and log g = 8.14, the latter T cff = 
11,600 K and logg = 8.05. Using the evolutionary models 



of Wood (1994) with thick surface H layers (~ 10 M*) our 
values translate into a mass of 0.64 Mq. 

As impressive as the fit in Figure ft] is, the discrepancies 
between model and data are still dominated by systematic 
effects, rather than stochastic noise. This makes the values 
of x 2 we calculate useless for evaluating the error in our 
temperature and gravity determination; the error is domi- 
nated by real differences between the models and the mea- 
surements. These differences might arise from a variety of 
sources. 

One possibility is that the normalisation of our data 
to the continuum points was affected by errors in our cal- 
ibration of the instrumental response. Another possibility 
is that errors arise from the presence of metal lines not in- 
cluded in the models. The shape of H8 = Hc( and He are 
probably affected by depression of the intervening contin- 
uum by the Ca II A3933 line and He is contaminated by 
Ca II A3968. However, the metal lines cannot entirely ac- 
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Figure 2. Location of best fits to the individual Balmer lines in 
the log g-T e fj plane. Fits to individual lines are denoted by the 
line designations; the best fit using all the lines (see Figure |lj) is 
shown as a filled circle. 



will provide the opportunity to improve atmospheric mod- 
els, perhaps even to infer the correct temperature profile. 
Furthermore, individual spectra in our time series, which 
themselves have a signal to noise ratio of fOO, span a tem- 
perature range of about 500 K and will show how the at- 
mospheric structure should change with model temperature. 
This information about the derivative of the thermal profile 
should provide extremely valuable constraints on model at- 
mospheres. Our data are available upon request to anyone 
interested in exploring these problems. 

The final factor that might affect our model fit is the 
presence of relatively large amplitude pulsations. We have 
shown in Paper I that the velocity fields associated with the 
pulsations of G 29-38 are detectable in our data. These mo- 
tions alter the average spectrum, mainly by Doppler broad- 
ening the absorption lines very slightly. The models we will 
describe in §^ allow us to incorporate this effect into the 
model spectrum. Fitting the data with these velocity broad- 
ened spectra yields a slight improvement in the quality of 
the fit to the average spectrum, but an insubstantial change 
in the values of log g and T e ff we infer from the models. 



count for the problem, because the fits to H/3, H/y, and H<5, 
which should be unaffected by metals, are also less than 
perfect. 

The discrepancy in the fits might also be d ue to a prob- 
lem discussed by Koester, Allard, & Vauclair (L994). They 



3 LINE PROFILE VARIATIONS 



found that even the best prescription for convection yields 
a temperature structure in model atmospheres that is only 
approximately correct. Consequently, the synthetic spectra 
the models produce will not match observed spectra at ev- 
ery wavelength simultaneously nor will they match all the 
line profiles, which are highly sensitive to the run of tem- 
perature with depth. To explore this possibility, and to see 
how our final fit is affected by individual lines, we have fitted 
separate models to each of the Balmer lines. 

In Figure ^, we plot the values of log g and T ett - for the 
best-fitting models of the individual lines H/3 through H9. 
They span a range of temperature from If ,620 to f2,885 K. 
With the exception of H9 = H?7, the fits show a trend to 
higher T e g and lower logp as the excitation level increases. 
The location of the fit using all lines compared to the in- 
dividual lines shows that the gravity fit is dominated by 
the lines of higher excitation, which are known to be more 
gravity sensitive. Conversely, the temperature is fixed by the 
lower excitation lines, which change more rapidly with T e g. 

Figure |^ is not meant to suggest that the individual 
fits give answers inconsistent with the global fit, rather they 
should be regarded as helping to establish the uncertainty of 
the temperature and gravity determination. It is clear that 
high signal-to-noise ratio alone is not enough to improve the 
determination of temperature and gravity. Better calibration 
of the instrument or improvements to the models (or both) 
will be necessary. The systematic trend in the individual 
Balmer line fits suggests the problem lies with the models, 
as do the results of the following section. There we compare 
fractional amplitudes calculated from data and models, and 
find further discrepancies. In the fractional amplitudes, cal- 
ibration errors should cancel to first order. 

If the models turn out to be responsible for the discrep- 
ancies we measure, our high signal to noise average spectrum 



Robinson et al. (1995) have described and implemented a 
method for distinguishing the value of £ in ZZ Ceti stars 
by comparing pulsation amplitudes at different wavelengths. 
Their method relies upon the increased importance of limb 
darkening in the ultraviolet. Non-radial pulsation modes of 
every £ suffer from geometric dilution in their amplitudes 
when averaged over the visible hemisphere. This dilution in- 
creases with I as cancellation between regions with opposing 
phase becomes important. However, at short wavelengths 
increased limb darkening diminishes the effect of the cancel- 
lations. For £ < 3, the net result is that mode amplitudes 
increase in the ultraviolet relative to thei r opt ical values in 
a way that depends on £. Robinson et al. (1995) used optical 



and ultraviolet high speed photometry to measure this ef- 
fect for the 215 s pulsation mode in G I I7-BI5A. They were 
able to conclude that this mode is £ = I. Fontaine et al. 



( 1997 ) reanalyzed the same data using independent models 
and arrived at the same conclusion, although they differed 
with Robinson et al. in the model temperature that best fits 
the data. 

The models Robinson et al. ( 1995 ) calculated also 
showed ^-dependent differences in the behaviour of pulsa- 
tion amplitudes within the absorption lines. These too arise 
from the effects of limb darkening on modes of different £. 
Figure ^ shows these changes for low values of £, which are 
expected to dominate the modes observed in white dwarfs. 
We calculated these curves using a modified version of code 
ori ginall y provided to us by E. L. Robinson (see Robinson et 
al. 1995). Instead of integrating over a limb darkening law, 



we integrate over intensities tabulated for different values of 
fi, as described in Consequently, Equations 3a and 3b 
in Robinson et al., which represent the equilibrium flux and 
the flux changes due to the pulsations, become: 



F x = 2kB£ I h(g,T ,fi)tidti, 



(1) 
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Figure 3. Models for the wavelength dependent flux varia tions 
for modes of t = 1 through 4, after Robinson et al. ( 1995| ), We 
have used the values of log g and T c g from our best fit to the av- 
erage spectrum and convolved the output with a Gaussian of 5 A 
width to match our Keck observations. All curves are normalised 
at 5500 A. The solid line is i = 1, the short dashed line is i = 2, 
the dotted line is I = 3, and the long dashed line is I = 4 
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Here, Ro is the equilibrium stellar radius, To the equilib- 
rium temperature, e is the amplitude of the radius changes 
induced by the pulsations, and |^ is the Lagrangian deriva- 
tive of temperature with respect to radius. Together Pe(fl) 
and kim represent the surface distribution of the tempera- 
ture changes after integrating in the (f> direction. Pe(fJ>) is a 
Legendre polynomial depending only on fj,, and kt m depends 
on the angle between the pulsation axis and the observer's 
line of sight . Our notation is slightly different from Robinson 
et al. (1995) in that we have separated the time dependence, 
cos (at), from ke m . Finally, while I\(g,T,[i) comes directly 
from the tabulated models, 91 x must be calculated by 

taking differences be tween models of different T e g. 

Robinson et al. (1995) have emphasized the useful prop- 
erties of the ratio AFx /Fx for I identification. It is not sen- 
sitive to the flux calibration of the data and varies with 
wavelength in a way that does not depend on mode inclina- 
tion and m (see also the Appendix). In Figure ^ we show 
AFx /Fx for modes of I = 1 through 4. The curves have been 
normalised to one at 5500 A. We have not included i = 0, 
which resembles 1=1 with slightly smaller modulations, 
because it is clear from the long periods of the modes that 
they cannot be radial pulsations. 

In order to compare our data to the models in Fig- 
ure ^, we have fitted the amplitudes and phases of the 11 
largest pulsation frequencies in each 2 A wavelength bin us- 
ing the function Acos(2-7r ft — (f>). During these fits, we held 
the mode frequency, /, constant at the values tabulated in 
Paper I, and fit the 11 amplitudes and phases simultane- 
ously. Of these 11 modes, five are combination frequencies, 
i.e., frequencies which have values that are sums or differ- 
ences of larger modes. We have discussed the nature of com- 



bination frequencies in Paper I and will return to them at 
the end of this section. Figure ^ shows the fractional ampli- 
tudes at each wavelength for the six physical modes and for 
the four largest combination frequencies. The fifth, F3— Fl 
at 10,322 s, we judged too noisy to include. The qualitative 
similarity between the data and the models is striking. Even 
in the line cores the models predict the behaviour of mode 
amplitudes quite well, despite the fact that within ~ 1 A of 
the core non-LTE effects in the atmosphere are important, 
which are not included in the model. We emphasize that 
these are predictions in the literal sense; Koester calculated 
the atmosphere models before we acquired the data. 

In the models we have discussed so far, the pulsations 
have the same phase at every wavelength. However, the 
phases of the physical modes, shown in Figure |^, show dis- 
tinct changes in the vicinity of spectral features. For mode 
Fl, the phase changes within absorption lines bear the sig- 
nature of a velocity induced variation; they change with the 
derivative of the spectrum. We will return to these phase 
changes in when we have models capable of reproduc- 
ing them. In addition to the phase changes within the lines, 
Figure |^ also shows a small slope in the continuum phases, 
indicating that pulse maximum in blue light arrives earlier 
than in red by a few seconds. This slope is not reproduced 
by our models, and remains a mystery. 

Careful inspection of the amplitudes of the real modes 
in Figure ^ reveals that those for mode F4, at 776 s, show a 
different shape than for the other modes; they increase more 
sharply in the line cores and curve more steeply in the con- 
tinuum. F4 was already noticeably different from the other 
real modes in Paper I, where we found it had a larger ve- 
locity to light amplitude ratio than any of the other modes, 
and it produced a stronger harmonic. 

To help see the differences in F4 we have normalised 
the amplitudes of Fl and F4 at 5500 A and plotted the 
modes together in Figure ^. The most likely explanation for 
the differences we see is that F4 has a value of I different 
from Fl. Qualitative application of the theoretical models 
demands that £ must be higher for F4, to yield the larger 
changes in amplitude we observe. Furthermore, the large 
contrast between modes of £ = 2 and higher in the models, 
versus the more modest differences seen in Figure [], point 
toward the conclusion that F4 is £ — 2 and Fl is I = 1. 

To test this possibility further, we can make a direct 
comparison between the data and models. In Figure ^ we 
have plotted the theoretical and observed amplitudes for 
modes 1, 2, 4, and 5. To establish the normalisation, we 
multiplied the theoretical curves by the amplitudes of each 
mode at 5500 A. We have used a model with gravity and 
equilibrium temperature inferred from our fit to the aver- 
age spectrum in §[]. As with the average spectrum, Figure ^ 
shows discrepancies between the data and the models. The 
slope of the amplitude changes is steeper in the data than in 
the models, worsening the fit at short wavelengths. Correct- 
ing this would require hotter models at shorter wavelengths, 
the same trend required to fit to individual Balmer lines in 
the average spectrum. We could not find a model at any sin- 
gle temperature that offered a substantially better fit than 
the model we have used. 

In spite of the difficulties with the fits, it is clear that 
£ — 1 is a better match to modes Fl, F2, and F5 while I = 2 
is a better fit to mode F4. Interestingly, this is most appar- 
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Figure 4. Wavelength dependent amplitudes for the 6 largest 
modes and 4 largest combination frequencies in G 29-38. To ease 
comparison, the same logarithmic scale was used for all panels. 
Note the small peak seen in the amplitudes of Fl at 3933 A, 
corresponding to the Ca II line. 



ent in the continuum variations between 4500 and 4700 A; 
the models are too poor a fit within the lines to provide 
a measure of £ there. Nonetheless, our original expectation 
that the line profiles would offer the most sensitive £ dis- 
criminant are borne out by the large differences between 
the amplitudes at line center in F4 and those in the other 
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Figure 5. Wavelength dependent phases for the 6 largest normal 
modes in G 29-38. 
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Figure 6. Amplitude comparison for Fl (solid line) and F4 
(discrete points). The amplitudes have been normalised to 1 at 
5500 A. 



modes. We have included F5 as the representative of the 
noisier low amplitude modes. Even with the higher noise it 
is clear that £ — 1 is a better fit to this mode, as it is to 
every mode except mode F4. 

Even though the quantitative agreement is poor at 
many wavelengths, the presence of two apparently different 
values of £ makes our identification secure. The values of £ 
chosen are independent of the model temperature. Lowering 
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Figure 7. Comparison to models for modes 1, 2, 4, and 5. The 
solid line is i = 1, the dashed line is I = 2. The models were 
calculated using logg = 8.05 and equilibrium T e g = 11850 K, as 
derived from the fits to the average spectrum. 



the effective temperature of the model chosen diminishes the 
changes in mode amplitudes, so that at some T, the £ — 2 
model would fit mode Fl and those like it, but at that tem- 
perature, no value of £ fits F4. Likewise, attempts to fit F4 
with £ — 1 by increasing model temperature leave no I of 
lower value to fit the other modes. Consequently, on the ba- 
sis of qualitative behaviour alone, we can conclude that the 
776 s mode is £ = 2, and the remaining 5 modes are l = \. 

Finally, we consider the combination modes shown in 
Figure [l|. The combination frequencies we see in G 29-38 
and other ZZ Ceti stars probably do not arise from eigen- 
mode pulsations, but from a non-linear transformatio n of the 



modes in the outer layers of the white dwar f (Brickhill 1992b 
cf. Brassard, Fontaine, & Wesemael 1995). This mixes the 
modes present, generating signals at sums and differences 
of the mode frequencies. The £ character of these modes 
depends on the £ of the modes which produce them. For 
example, we expect the combination of two (£,m) = (1,0) 
modes to have £ — and £ = 2 components, while those of 
two (1,-1) modes or two (1,1) modes should produce only 
I = 2. These expectations arise from the mathematical prop- 
erties of spherical harmonics only. It is impossible to make 
quantitative predictions about combination modes without a 
detailed theory explaining how they are produced. We have 
shown two combination frequencies along with models in 
Figure [| Like the other combinations, they most resemble 
the modes we have identified as £ — 1. We have already 
discussed the implications of this in Paper I. 
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Figure 8. Comparison to models for frequencies F1+F2 and 
F2+F2. The solid line is I = 1, the dashed line is I = 2. 



4 MODELS INCLUDING VELOCITIES 

Apart from the problems with quantitative fits to the model 
atmospheres, there is a significant qualitative difference be- 
tween the observed and model amplitudes in Figure [jj; the 
observed amplitudes are asymmetric within the absorption 
lines, while the model amplitudes are not. We showed in Pa- 
per I that the Balmer lines not only show changes in flux 
during a pulsation cycle, but also changes in line-of-sight 
velocity. In that paper, we treated these as separable com- 
ponents of the spectral variations, but in reality they are 
components of a more complex line profile variation. In an- 
ticipation of improvements in the model fits, in this section 
we will calculate these variations by incorporating veloci- 
ties into the flux integrals we calculate at each wavelength. 
This requires that we follow a mor e general treatment than 

to produce Equation |l[ 
Nather 



that used by Robinson et al. ( 199 



Our d evelopment will rely upon Robinson, Kepler, 
and use the same notation where possible. 



The effect of the velocity field on our calculations of 
the integrated flux is that the values of 7a we look up from 
the tabulated spectra, and the values of we calculate 
from them, must be chosen using wavelengths adjusted for 
the velocity of each surface element. Thus I\ in Equation |l|, 
and in Equation [| are now functions of velocity. Since 
the velocity has a time dependence, both quantities are now 
also implicit functions of time. We may write the total flux 
at each observed wavelength as a function of time as: 



F(X ohs ,t) = F To (A obs ) + AF T (X b S ,t), 



(3) 



where the Tb is a reminder that the first term comes from 
the model at equilibrium temperature, and the T indicates 
that the flux changes in the second term are due to changes 
in tempe rature , which dominate all other sources (Robin- 
son et al. 1982). It will help in understanding Equation ^ to 
think about two different limits. The first is a (hypothetical) 
pulsation mode which has surface motions but no tempera- 
ture changes. For these modes, the second term would be 0. 
The first term would be a function of time for any A b s near 
spectral features that could be Doppler shifted in and out of 
Aobs by the changing velocities. The second limit to think of 
is a mode with temperature but no velocity variations, then 
the two terms in Equation ^ reduce to the expressions given 
by Equations [j] and [| An important feature of Equation ^ 
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that is not present in Equations [I] and |^ is the possibility 
for the phase of flux maximum to differ from the phase of 
temperature maximum. This means the time of flux maxi- 
mum can differ with wavelength, a possibility that was never 
allowed by Equation |l|. 

Our expressions for the terms in Equation ^ must also 
be more general than before. We can write Ft generally as 



Ft (A bs) — Rq 



I(g,T ,iJ,,\) fj,d^< 



(4) 



which is analogous to the expression used by Kepler (1984) 
in his discussion of line profile variations due to r-mode pul- 
sations. As we have noted, the time dependent velocities 
enter into the wavelength A, so -Pr (A o bs) is a function of 
time. Likewise, AFt is given by 



AFr(A b s , t) 



(R ^)eR 2 e-^ 
or 



1 dI(g,T,^,\) 



OT 



T 



(•>) 



where £ r is the assumed functional form for the perturba- 
tions in stellar radius, 



(6) 



Yem(<d,&) is the spherical harmonic of degree I and order 
m in the coordinate system (0,$) aligned with the pulsa- 
tion axis, and a is the pulsation frequency. Compared with 
Equation || we have reverted to an expression in which the 
imaginary parts of the temporal and spatial dependenc es are 
included; moreover, we have followed Robinson et al. (f982) 
in adopting an extra term e~'^ to allow for the nonadiabatic 
effects which may introduce a phase difference between the 
radial displacement and the flux changes^]. The negative sign 
in the exponent indicates that for positive values of ip, max- 
imum flux lags maximum radial displacement. 

In Equations [| and ^, the value of A should be the wave- 
length from which light is Doppler shifted into A b s , or, 



A = A n 



Vrad(Af><M) 



(7) 



to first order. We use u ra d to express the velocity component 
projected into our line of sight. Calculating f ra d requires 
first an expression for the pulsation velocities at the stellar 
surface. 

In the frame of reference (0,$) aligned with the pul- 
sation axis, the compon ents o f the pulsation velocities are, 
following Dziembowski (1977), 



V r = 
V B = 

V* = 



ie\g\ dY em (Q,$) 



a 



and 



a sin 



(8) 



In Robinson et al. (1982), Equations 20, 23 and 24 should not 



have negative signs on their right hand sides. The quantity 44- is 
positive for adiabatic pulsations, since maximum radial displace- 
ment corresponds to maximum temperature. This has no effect 
on any of the conclusions Robinson et al. presented. 



For g-mode pulsations, V r is small compared to the other 
two and can be ignored. 

To get u ra d requires transforming the remaining veloc- 
ity expressions into our reference frame and projecting them 
along our line of sight. Then Equations ^| through (7] are all 
we need in principle to calculate colour-dependent pulsation 
amplitudes in the presence of non-zero velocities. In prac- 
tice, this would be cumbersome and inefficient for arbitrary 
inclination and m, so we have further simplified the prob- 
lem by aligning the pulsation axis with our line-of-sight and 
holding m — 0. Contrary to initial expectations, this sim- 
plification comes at almost no expense; the results for this 
case scale easily to arbitrary choices for inclination and m. 
In the Appendix we demonstrate that this is true and show 
how the scaling is done. 

This simplification allows us to dispense with the ex- 
pression for Vi> and to replace and $ in all equations 
with 9 and <j> for our reference frame. The <j> dependences of 
^rad, £r, and I\ also disappear, and the spherical harmonics 
reduce to Legendre polynomials. Finally, to make our fits 
independent of the particular choices for stellar radius and 
pulsation frequency, we have introduced the fitting parame- 
ters: 



ax 



Ro 
£9 



5T 
6r 



and 



(9) 



Thus the final expressions we evaluate numerically are: 



Ft (A bs) 
AF T (A obs ,i) 



27ri?o 



I{9,T ,ij,, A oba [l 



«rad 



) A* dn, 



2-kRq ar cos(at — tp) x 
1 dI{g,T,f,,X oh3 [l-^] 







Urad 



(M) 



dPi(n) 
1 d6 



OT 



) sin(o-t). 



Pt(fj.) ndn, 



T 



(10) 



Where we have kept only the real parts of the temporal 
variations. The free parameters in these equations are the 
temperature and gravity of the equilibrium model, To and g, 
the (non-adiabatic) phase shift ip, and the amplitudes given 
by Equation ^. The equilibrium radius, Ro, cancels when we 
calculate the fractional amplitude. 

We have modified the code originally provided by E. 
L. Robinson to perform these integrals at a series of time 
steps covering one pulsation cycle. Then we have convolved 
the output spectra with a Gaussian to emulate seeing. Fi- 
nally, we have calculated the pulsation amplitude, phase, 
and mean spectrum at each wavelength. 

In Paper I, we defined the quantity A$v as the phase 
difference between maximum light and maximum velocity. 
In the formalism of this paper, 



A<Ev = ip- 



(11) 



For the ip = case, where flux and radial displacement are 
in phase, the ^ delay enters because of the time derivative 
used to get u ra d- Positive values of ip then reduce A$y- by 
delaying the maximum light so that it arrives less than -| be- 
fore maximum velocity. In all of the modes for which we can 
detect velocities, A$y lies in the first quadrant, implying 
that flux maximum is delayed compared to the adiabatic 
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Figure 9. Fits to Fl and F4 including velocities. The amplitudes 
listed are those described by Equation ^ in the text. 



case (or velocity maximum advanced, which is harder to 
imagine). This is a profound result, and the first direct ob- 
servational constraint on the behaviour of eigenmodes near 
the surface of a white dwarf star. 

Comparing our models to the data presents something 
of a challenge, because of the generally poor match for any 
choice of amplitude and phase. We have simplified the prob- 
lem by choosing the equilibrium temperature and gravity 
from our fit to the average spectrum. Then we have fit- 
ted (it directly, by insisting that the fractional amplitude 
of the flux changes in our model match those in the data 
at 5500 A. To choose appropriate values for a v and ip, we 
used our model to calculate several time series of synthetic 
spectra and then we reduced them in the same way as the 
real data. This allowed us to calibrate the scaling introduced 
by the integrals in Equation [l(] for various choices of £, a v , 
and ip. Then we used these scale factors to generate models 
based on the properties measured for modes in Paper I. The 
size of the velocities we use to match our data i s quite similar 
to the size predicted by Robinson et al. (1982), ~ 7 kms -1 . 

Figure ^| shows the results of these calculations for the 
H/3 line in modes Fl and F4. Our models now reproduce the 
asymmetry in the amplitude plots, and show phase changes 
within the line profiles. These phase changes could never 
be reproduced without including velocities in the models. 
However, the overall quality of the amplitude fits within 
lines is not much improved, and is not changed at all in the 
continuum. The phases are a good fit within spectral lines 
for the larger modes, but our models still do not reproduce 
the slow change in phase observed in the continuum (see 
Figure ^ . For now, this remains a mystery, but probably a 
mystery hiding interesting pulsation physics. 



Finally, we consider the behaviour of the ratio between 
the observed velocity and flux amplitudes. As we show in the 
appendix, this ratio does not depend on the inclination or 
m of the modes, and therefore cannot be used to constrain 
those quantities. However, it is sensitive to the value of £ as 
evident from the different forms for flux and velocity in the 
integrals of Equation [l^. The flux depends on the Legendre 
polynomial while the velocity depends on its derivative and 
in addition is weighted towards the limb by the projection 
onto the line of sight. According to our models, for a fixed 
value of — , an I = 2 mode should have ~ 4 times higher 
apparent velocity to light ratio than an £ — 1 mode. It is 
then no surprise that mode F4, which we have identified 
as £ = 2, has the highest velocity to light amplitude ratio 
(see Paper I) of any of the modes we have detected, a result 
which supports our identification of that mode as £ — 2. 



5 SUMMARY AND CONCLUSIONS 

Asteroseismology of ZZ Ceti stars has been impeded for al- 
most two decades by the lack of a reliable method for mode 
identification. We have tested a new method that uses high 
signal to noise time-resolved spectroscopy to measure the 
wavelength dependence of optical pulsation amplitudes. We 
have found that the ^-dependent changes in amplitude pre- 
dicted by the models are also present in our data, allowing 
us to assign values of £ to six modes in the star G 29-38. 

Initially, this will be of greatest benefit to seismologi- 
cal models of G 29-38. None of the modes we have identi- 
fied is short enough to yield immediate constraints on the 
mass of t he H layer, as is possible for G 117- B15A (Robin- 
son et al. |l995|) , GD 1 65 (Bergeron et al. |l993| ) , and G 226-29 
(Fontaine et al. |l992j ). However, a concerted effort to match 
the six modes of known £ may yield a unique solution. If 
not, it is possible to return to this star with the hope that 
four more hours of data will allow identification of a differ- 
ent set of modes. Sooner or later, we will have a definitive 
asteroseismological solution for this star; the main obstacle 
has been removed. 

Our method can also be extended to other ZZ Ceti stars, 
although the fainter, low amplitude stars will require longer 
runs. We have already begun a programme to identify £ in as 
many stars as is practical in the observing time available to 
us. This programme should allow us to measure structural 
properties of enough stars to answer some long-standing 
questions about the DA stars, such as the masses of the H 
and He surface layers, whic h have been the subject of so me 
controversy (Shipman 1996; Fontaine & Wesemael 1996). 



Apart from the impact our results will have on ZZ Ceti 
seismology, they have also opened new windows into the sur- 
face physics of pulsating white dwarfs. We have measured 
two new diagnostics of the behaviour of pulsations near the 
photosphere: the amplitude of the pulsation velocities, and 
the phase lag between flux maximum and velocity maxi- 
mum. These can be compared quantitatively with predic- 
tions of nonadiabatic p ulsati on theories (e.g., Lee & Bradley 
199!: ; Wu & Goldreich 199S ) , as well as with measurements 
for other stars. We have already detected velocity changes 
in one other star (HL Tau 76). 

Finally, our high signal to noise measurements provide 
a serious challenge to the atmospheric model fits. There are 
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significant discrepancies between the models and our average 
spectrum and between the models and the spectral changes 
caused by pulsation. We have not identified the source of 
these differences, but hopefully our data and future data 
like them will provide assistance to modelers in their con- 
tinuing attempts to understand the atmospheres of white 
dwarf stars. 
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APPENDIX A: ARBITRARY INCLINATION 
AND VALUE OF m 

The usual way to transform an arbitrary spherical harmonic 
into the coordinate system aligned with our line of sight is to 
recognise that the 2£+l spherical harmonics of a given I form 
a complete basis set, so it is always possible to express an 
arbitrary spherical harmonic as a sum of spherical harmonics 
in our coordinate system, 



E 



R n 



zY lm ,(9, 0), 



(Al) 



The coefficients R„ 



1977 



are a function of the inclin ation angle 
Robinson et al. [1982J). 



Then the 



(see Dziembowski 

m' 7^ terms of the sum cancel to zero in any integration 
over the visible hemisphere, so that only the m! = term 
aligned with our line of sight remains. 

This will not work directly for the expressions we inte- 
grate in Equations ^ and 0, because they contain I(g,T, fj,, A) 
and its derivative with respect to T, which may be arbitrar- 
ily complex functions of the eigenmode velocities. However, 
if we expand 7 in a Taylor series about A, 



+ 



a/( 5 ,T, M ,A) 



dX 



(A — A b s ), 



(A2) 



do the same for ^ , and use Equation |^, our expression for 
the time dependent flux becomes 



F(X obs ,t) 



+ 



2ttRI 



D 2 A b s 
H 



K a,Te e 



/(m> <^obs) H dfj, 



r)X 



/i d/i 



dI(n,X ohB 



dT 



H dfi. (A3) 



r 



We have suppressed the g and T dependencies, which van- 
ish upon choosing an equilibrium model. We have also left 



out the cross term which includes the product of J|r and 
|j because these are both small quantities. The remain- 
ing expression contains separate terms for the equilibrium 
flux (Fo), the flux changes due to velocity shifts (AF V ), and 
the flux chan ges due to temperature changes(AFr). Using 
we can express both « ra d and Yi m in the above 



Al 



Equation 

integral into the coordinate system (0, <f>); for both, the terms 
with m ' integrate to zero in the <f> direction (Dziem- 
bowski 1977), so only the m — term aligned to our line of 



sight remains. Thus we may write: 

F(X bsyt) = i^O (A bs) 

+ Rom&F v (X obs ,t,i = 0,m = 0) 
+ Ro m AF T (X ohs ,t,i = 0,m = 0). 



(A4) 
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This expression permits us to relate AF/F(\ \, B , i 7^ 
0,m =f 0) to AF/F(Aobs,i = 0, m = 0) via a simple scal- 
ing factor i?om- As Ro m does not depend on A b s , the results 
shown in Figure 3 are general for arbitrary inclination and m 
value. Also, because the factor is the same for both the tem- 
perature and velocity induced flux changes, the observed ra- 



tio a v /a,T does not depend on i and m. Finally, Equation A.4 



shows that the line-of-sight velocities and flux variations arc 
separable to first order, a result used implicitly in Paper I. 

It is important to recognise where these approximations 
may break down. When the second derivative of the spec- 
trum dominates the first, which can happen in the central 
few Angstroms of an absorption line, our expansion to first 
order in A is not reliable. 
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